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ABSTRACT 

We report the discovery, from WASP and CORALIE, of a transiting exoplanet in a 1.54-d orbit. 
The host star, WASP-36, is a magnitude V = 12.7, metal-poor G2 dwarf (Teff = 5959 ± 134 K), with 
[Fe/H] — —0.26 ±0.10. We determine the planet to have mass and radius respectively 2.30 ±0.07 and 
1.28 ± 0.03 times that of Jupiter. 

We have eight partial or complete transit light curves, from four different observatories, which 
allows us to investigate the potential effects on the fitted system parameters of using only a single 
light curve. We find that the solutions obtained by analysing each of these light curves independently 
are consistent with our global fit to all the data, despite the apparent presence of correlated noise in 
at least two of the light curves. 

Subject headings: planetary systems - planets and satellites: detection ~ planets and satellites: fun- 
damental parameters - stars: individual (WASP-36) - techniques: photometric 



1. INTRODUCTION 

Of the 171 confirmed transiting planetary system^, 
the majority have been discovered from the g round, from 
surv eys such as WASP (|Pollacco et all 120061 ) and HAT- 
net (jBakos et al.ll2004) . Although the Kepler space mis- 
sion is discovering an increasing nu mber of planets and 
even more can didate planets (e.g. iBorucki et al.l [20Tol : 
iBorucki et al.l [2011.). the ground-based discoveries have 
the advantage that the host stars are generally brighter. 
This allows radial-velocity measurements to measure the 
planetary mass, and is conducive to further characterisa- 
tion observations, such as measuring occultations in the 
infrared to probe atmospheric temperature and struc- 
ture. 

Many of the current questions in exoplanet science are 
being addressed by analysing the statistical properties 
of the growing ensemble of well characterised transiting 
planetary systems. Here we report the discovery of a 
transiting planet orbiting the V ^ 12.7 star WASP-36 {— 
2MASS J08461929-0801370) in the constellation Hydra. 

2. OBSERVATIONS 
2.1. WASP photometry 
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WASP-36 was observed in 2009 and 2010 by WASP- 
South, which is located at the South African Astro- 
nomical Observatory (SAAO), near Sutherland in South 
Africa, and by SuperWASP at the Observatorio del 
Roque de los Muchachos on La Palma, Spain. Each in- 
strument consists of eight Canon 200mm f/1.8 lenses, 
each equipped with an Andor 2048 x 2048 e2v CCD 
camera, on a single robotic mount. Further details 
of the instrument, survey and data reduction proce- 
dures are described in .Pollacco et al] (|2006f ) and details 
of the candidate selection p roce dure can be found i n 
iCollier Cameron et al.l (|2007l) and lPollacco et al.l (|2008D . 
A total of 13781 measurements of WASP-36 were made 
between 2009 January 14 and 2010 April 21. 

WASP-South 2009 data revealed the presence of a 
transit-like signal with a period of '-^ 1.5 days and a depth 
of ~ 15 mmag. The WASP light curve is shown folded 
on the best-fitting orbital period in Figure [1] 

2.2. Spectroscopy 

Spectroscopic observations of WASP-36 were made 
with the CORALIE spectrograph of the 1.2-m Euler- 
Swiss telescope. Simultaneous spectra of a thorium- 
argon emission line lamp were obtained in order to cal- 
ibrate the stellar spectra. A total of nineteen spectra 
were taken between 2010 March 11 and 2011 January 
11, and processed using the standard CORALIE data 
reduction pipeline (jBaranne et al.llT996() . The resulting 
radial- velocity data are given in Table [Tl and plotted in 
Figure [5] In order to rule out non-planetary causes for 
the radial-velocity variation, such as a blended eclips- 
ing binary system , we examined the bisector spans (e.g. 
iQueloz et al.ll200"ll ). which exhibit no correlation with ra- 
dial velocity (Figure. [2]) , as expected. 

2.3. Follow-up photometry 

We have a total of eight high-precision follow-up light 
curves of the transit of WASP-36b, summarised in Table 
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TABLE 1 








Radial- 


VELOCITY (RV) AND LINE 


BISECTOR SPAN 


(BS) 






MEASUREMENTS OF WASP-36 




BJD(UTC) 


RV 






BS 


-2 450 000 


(km s^i) 


(km s'^) 


(k 


m s"-*-) 


5266.6926 


— 12.843 


0.027 




0.079 


5293.5807 


-13.591 


0.021 




-0.017 


5304.6409 


-13.212 


0.057 




0.025 


5305.5493 


-13.423 


0.025 




-0.046 


5306.6464 


— 12.847 


0.025 




0.014 


5315.5600 


-12.981 


0.023 




0.013 


5316.5339 


-13.624 


0.028 




-0.005 


5317.5642 


-12.977 


0.033 




-0.046 


5320.4833 


— 12.821 


0.030 




-0.034 


5359.4568 


-13.527 


0.040 




0.002 


5547.8347 


-12.842 


0.024 




0.007 


5561.8393 


-12.866 


0.026 




0.003 


5562.8651 


-13.306 


0.025 




-0.018 


5563.8184 


-13.402 


0.033 




-0.070 


5564.7341 


-12.920 


0.024 




-0.011 


5565.8211 


-13.467 


0.029 




0.056 


5567.8088 


-12.870 


0.023 




0.001 


5570.7984 


-12.931 


0.023 




-0.028 


5572.7427 


-12.944 


0.023 




-0.068 



O Differential aperture photometry was performed us- 
ing the IRAF/DAOPHOT package for T RAPPIST and 
FTN data, and the ULT RACAM pipehne HDhillon et alj 
I2007t IBarros et aLll201l'l ) for the LT data, with aperture 
radii and choice of comparison stars optimised to give 
the lowest RMS of the out-of-transit photometry. 

3. DETERMINATION OF SYSTEM PARAMETERS 
3.1. Stellar parameters 

The individual CORALIE spectra of WASP-36 were 
co-added to produce a single spectrum with a typ- 
ical S/N of around 60:1. The standard CORALIE 
pipeline reduction products were used in the analy- 
sis. The sp ectral analysis was perform ed using uclsyn 
(|Smith fc D woretskv 1988; Smithlll992[) and the methods 
given in [Gillon et al.. (■200ft) . The parameters obtained 
from the analysis are listed in Table [3l The elemen- 
tal abundances were determined from equivalent width 
measurements of several clean and unblended lin e s. Th e 
lines used are t hose l isted i n iGonzalez fc Laws! (I2000D . 
iGonzalez et"all (|2001[ ). and iSantos et al.l (|2004alV A 
value for microturbul ence, £t, was d etermined from Fe i 
using the method of iMagainl ()1984l ). The quoted error 
estimates account for the uncertainties in Tea, logg and 
^t, as well as for the scatter due to measurement and 
atomic data uncertainties. 

The projected stellar rotation velocity, wsini, was de- 
termined by fitting the profiles of several unblended Fe i 
lines in the wavelength range 6000-6200A, using the ro- 
tation broadening function of (jGravl [20081 Ch. 18). We 
used an instrumental FWHM of 0.11 ± 0.01 A, deter- 
mined from the telluric lines around 630 nm. The mea- 
sured V sini is sensitive to the adopted value of Vmar. - The 
appropriate value of Umacis 4.0 kms~^ accor ding to ( Graivl 
12008', p. 507), but 2.9 kms'^ according to lBruntt eraD 
((2010,). These values imply vsini = 2.9 ± 1.3 kms""'^ 
and usini = 3.7 ± 1.1 kms^^ respectively. We take the 
weighted average of these two values as the best-fitting 
one, vsini = 3.3 ± 1.2 kms~^ . The quantity measured 



is approximately the quadratic sum of v sin i and Wmac 
4.9 ± 0.8 kms-i ). 

3.2. Neighbouring objects 

The Two Micron All Sky Survey catalogue 
(jSkrutskie et al.l |2006[) reveals the presence of four 
fainter stars close on the sky to WASP-36. There is 
no evidence from analysis of catalogue proper motions 
that any of these stars are physically associated with 
WASP-36. The stars are separated from WASP-36 by 
4", 9", 13"and 17", meaning that they fah well within 
the WASP photometric aperture, which has a radius of 
48"(3.5 pixels), but outside of the 1" CORALIE fibre. 

In the absence of reliable optical catalogue magnitudes 
for all of these objects, it was necessary to measure their 
fiuxes to quantify the effects of blending in the photom- 
etry. The fluxes were measured from images taken dur- 
ing the two transits observed with the 1.2-m Euler-Swiss 
Telescope (see Table [2]). The fluxes relative to that of 
WASP-36 are as follows: 0.012 (object at 4" separation 
from WASP-36), 0.00771 (9"), 0.00558 (13") and 0.00827 
(17"). Using these flux ratios, we corrected the WASP 
photometry to account for all four objects, and the high 
precision photometry to account for the object at 4", 
which is within the photometric apertures used. The 
magnitude of this correction is minimal, and had no sig- 
nificant (^ 1-cr) effect on the values of our best-fitting 
system parameters. 

3.3. Planetary system parameters 

CORALIE radial-velocity data were combined with all 
our photometry and analysed simultaneously using the 
Markov Chain Monte Carlo (MCMC) method. The best- 
fitting system parameters are taken to be the median 
values of the posterior probability distribution. Linear 
functions of time were fitted to each light curve at each 
step of the MCMC, to remove systematic trends. We 
use the current ve rsion of the MCMC code de s cribed 
in CoUier Camero n et al.l (|2007[ ): iPollacco et all (|2008l ) 
and Enoc h et al.l (|2010f l. The MCMC proposal parame- 
ters we use are: the epoch of mid-transit, Tc] the orbital 
period, P; the transit duration, T14; the fractional flux 
deficit that would be observed during transit in the ab- 
sence of stellar limb-darkening, AF; the transit impact 
parameter, b; the stellar reflex velocity semi-amplitude, 
Ki; the stellar effective temperature, T^s; the stellar 
metallicity, [Fe/H]; and ^/ecosw and ^/esinoj, where e is 
the o rbital eccentricity, and uj is the argument of perias- 
tron ([Anderson et al.ll2011[ ). The stellar mass was deter- 
mined as part of the MCMC analysis, using an em pirical 
fit t o [Fc/ H], Teff, and th e stellar density, (Enoch et al.l 
[20101: [Torres et al.l[2010[ ). 

The tr ansit light curves were m odelled using the formu- 
lation of iMandel fc Ago] ([2002D and limb- darkening was 
accounted for using a four-coefficient, non-linear model, 
employing coeffici ents appropria te to the passband from 
the tabulations of lClareti ([20001 12004D . The coefficients 
were determined using an initial interpolation in logg* 
and [Fe/H] (values from Table [3]), and an interpolation 
in T,^cff at each MCMC step. The coefficient values cor- 
respondi ng t o the best- fitting value of T^,,cB are given 
in Table 13.31 Because some of our photometry was ob- 
served in passbands not tabulated bv iClaret (,2000,. 2004) . 
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TABLE 2 

Observing log for follow-up transit photometry 



Light 


Date/UT 


Telescope / instrument 


Band 




icxp 


full or 


Seeing or 


Aperture 


Airmass 


curve 












partial 


defocus {") 


radius (") 


range 


(i) 


2010 Dec 13 


Euler^ / EulerCam 


Gunn r 


94 


120 


partial 


1.1 - 2.2 


4.3 


1.54 - 1.07 - 1.09 


(ii) 


2010 Dec 13 


TRAPPIST^/ TRAPPISTCAM 


clear 


756 


10 


partial 


3 


8.3 


1.84 - 1.11 


(iii) 


2010 Dec 25 


FTN^ / Spectral camera 


PS 2 


176 


60 


full 


4.3 


2.4 


1.52 - 1.14 - 1.22 


(iv) 


2011 Jan 02 


TRAPPIST / TRAPPISTCAM 


I + z 


296 


25 


full 


2 


6.4 


1.75 - 1.09 


(v) 


2011 Jan 05 


TRAPPIST / TRAPPISTCAM 


clear 


179 


18 


partial 


3 


7.7 


1.18 - 1.07 - 1.14 


(vi) 


2011 Jan 08 


TRAPPIST / TRAPPISTCAM 


clear 


269 


18 


partial 


3.5 


9.0 


1.10 - 1.44 


(vii) 


2011 Jan 15/16 


LT_^ / RISE_^ 


V + R 


1290 


9 


full 


6 


9.0 


1.95 - 1.25 


(viii) 


2011 Jan 21 


Euler / EulerCam 


Gunn r 


167 


60 


full 


0.45 - 1.0 


4.1 


1.46 - 1.20 



''1.2-m Euler-Swiss Telescope, La Silla, Chile 

'^Transiting Planets and Planetesimals Small Telescope, La Silla, Chile Ijjehin et al.ll20111 ; |http://www.astro.ulg.ac.be/Sci/Trappist| | 
'^Faulkcs Telescope North, Haleakala Observatory, Maui, Hawaii, USA 
''Liverpool Telescope, Observatorio del Roque de l os Muchachos, La Palma, Spain 
''Rapid Imaging Search for Exoplanets camera ({Steele et al.ll2003 : IGibson et al.ll200SI 1 



TABLE 3 

Stellar parameters and abundances from analysis of 

CORALIE SPECTRA. 



Parameter 


Value 


Parameter 




Value 




RA (J2000.0) 


08h46ml9.30s 


[Fe/H] 


-0 


26 


± 





10 


Dec {J2000.0) 


-08'= 01' 36.7" 


[Na/H] 


-0 


33 


± 





08 




5900 ± 150 K 


[Mg/H] 


-0 


08 


± 





08 


log9(cgs) 


4.5 ± 0.15 


[Si/H] 


-0 


17 


± 





06 




1.0 ± 0.2 kms-i 


[Ca/H] 


-0 


15 


± 





11 


V sin i 


3.3 ± 1.2 kms-i 


[Sc/H] 


-0 


11 


± 





12 


log A(Li;i^ 


1.69 ± 0.13 


[Ti/H] 


-0 


16 


± 





11 


Sp. Type 


G2 


[V/H] 


-0 


20 


± 





15 


Distance 


450 ± 120 pc 


[Cr/H] 


-0 


28 


± 





09 


Age 


1 - 5 Gy 


[Mn/H] 


-0 


44 


± 





10 


Mass 


1.01 ± 0.08 Mq 


[Co/H] 


-0 


19 


± 





12 


Radius 


0.94 ± 0.17 Rq 


[Ni/H] 


-0 


30 


± 





08 



Additional identifiers for WASP-36: 
USNO-Bl.O 0819-0221838 
2MASS J08461929-0801370 
ISWASP J084619. 30-080136. 7 



N ote: The sp ectral type was estimated from T^ff using the table 
of IIGravll200a. p . 507) . The mass and radius were estimated using 
the [Torres et al] 1)20101 ) calibration. 
''Microturbulent velocity 

''log A(Li) = log ^-^^^ + 12, where A^n and A^h are the number 
densities of Li and H respectively. 

we also tried using coefBcients corresponding to nearby 
passbands. None of our best- fitting system parameters 
was significantly affected by our choice of Claret pass- 
band; values changed by much less than their 1-a uncer- 
tainties. 

An initial MCMC fit for an eccentric orbit found e = 
0.012to oo8 = 55tii8 degrees), with a 3-cr upper-hmit 
to the eccentricity of 0.064, but we found this eccentric- 
ity is n ot signific ant. F ollowing the F-test approach of 
ILucv fc Sweeney! (jl97lD . we find that there is a 58 per 
cent probability that the apparent eccentricity could have 
arisen if the underlying orbit were actually circular. We 
therefore present here the model with a circular orbit, 
noting that the values of the other model parameters, 
and their associated uncertainties, are almost identical 
to those of the eccentric solution. 

We tried fitting for a linear trend in the RVs with the 
inclusion of an additional parameter in our MCMC fit. 
Suc h a trend (s uch as that found in the RVs of WASP- 
34, iSmallev et al.. ,2011.) would be indicative of a third 



TABLE 4 
Limb-darkening coefficients 

Claret band Light curves ai a2 as 04 

Cousins R WASP,i,ii,v,vi,viii 0.466 0.294 0.070 -0.128 

Sloan z' iii.iv 0.555 -0.099 0.348 -0.213 

Johnson V vii 0.389 0.523 -0.066 -0.082 



body in the system. The best-fitting radial acceleration 
is consistent with zero, indicating there is no evidence 
for an additional body in the system based on our RVs, 
which span 10 months. The orbital parameters we report 
are the result of a fit which does not allow for a linear 
trend in radial velocity. 

The system parameters derived from our best-fitting 
circular model are presented in Table[5] The correspond- 
ing transit and RV models are superimposed on our data 
in Figures [1] and [2l 

3.4. System age 

The measured v sin i of WASP-36 gives an upper limit 
to the rotational period, Prot — 14.4±5.9 days. This cor- 
responds to an upper limit on the age of ~ l.SlHGy 
using the gyrochronological relation of I Barnes! (|2007D , 
and a i?-magnitude of 13.3 derived fro m V=12.7 a nd 
B -V ^ 0.60 ± 0.04 (estimated using [Gravl [20081 p. 
507). 

In FigurelDw e plot WASP-36 alon gside the stellar evo- 
lution tracks of iMarigo et all ()2008[ ). From this we infer 
an age of 2.512 2 ^ly. The age determined from the 
lithium abunda nce of WASP-36 is poorl y constrained, 
but the work of iSestito fc RandichI (|2005f) suggests that 
the most likely age is 2 to 5 Gy. 

We searched the WASP photometry for periodic vari- 
ations indicative of star spots and stellar rotation, but 
no significant variation was detected. We place an up- 
per limit of 1.5 mmag at the 95 per cent confidence level 
on the amplitude of any sinusoidal variation. This null 
result is consistent with the low levels of stellar activity 
expected from a main-sequence G2 star. A lack of stel- 
lar activity is also indicated by the absence of calcium 
II H-l-K emission in the spectra. The uncertainties on 
the Call emission index, logi?^^(~ 4.5) are too large to 
allow meaningful constraints to be placed on the system 
ag e by using an activity - rotatio n relation such as that 
of lMamaiek fc Hillenbr'cffldl (|2008[ ). 

There is no evidence of any discrepancy between the 



4 



A. M. S. Smith et al. 



TABLE 5 
System parameters 



Parameter 



Orbital period 

Epoch of mid-transit 

Transit duration 

Ingress / egress duration 

Planet-to-star area ratio 

Transit impact parameter 

Orbital inclination angle 

Stellar orbital velocity semi-amplitude 

System velocity 

Orbital eccentricity (adopted) 

Orbital eccentricity (S-cr upper-limit) 

Stellar mass 

Stellar radius 

log (stellar surface gravity) 

Stellar density 

Stellar effective temperature 

Metallicity 

Planet mass 

Planet radius 

Planet density 

log (planet surface gravity) 

Scaled orbital major semi-axis 

Orbital major semi-axis 

Planet equilibrium temperature (uniform heat redistribution) 
System age (from figure |3]l 



The following constant values are used: AU = 1.49598 X 10^^ m; 



, = 1.89896 X 10 


2^ k, 


g; pj = 1240.67 kg 


; m 








TABLE 6 










Transit times 






Light curve 


E 


Tc 


O-Tc 


O-C 






(HJD, UTC) 


(min) 


(min) 


(i) 


-17 


2455543.70602 


5.72 


5.65 


(ii) 


-17 


2455543.70378 


1.18 


2.42 


(iii) 


-9 


2455556.00221 


0.62 


1.71 


(iv) 


-4 


2455563.68807 


0.33 


0.32 


(v) 


-2 


2455566.76718 


7.73 


6.63 


(vi) 





2455569.83686 


0.86 


-0.62 


(vii) 


5 


2455577.52412 


0.21 


-0.02 


(viii) 


9 


2455583.67344 


0.27 


-0.23 



ages derived from lithium abundance, gyrochronology 
and isochrone fitting. This suggests that the star has 
undergone httle or no tidal spin-up, despite the presence 
of a massive planet in a close orbit. 

3.5. Transit timing 

We measured the times of mid-transit for each of the 
eight follow-up light curves, by analysing each light curve 
separately, without any other photometry (see Section 
14. ip . The times are displayed in Table I3.3[ along with 
the differences, O — C between these times and those pre- 
dicted assuming a fixed epoch and period (Table [5]) . No 
significant departure from a fixed ephemeris is observed. 

4. DETAILED ANALAYSIS OF FOLLOW-UP LIGHT 
CURVES 

Because we have several follow-up light curves of 
WASP-36 from different telescopes / instruments, 
whereas many planet discovery papers rely on only a 
single such light curve, we take the opportunity here to 



Symbol 


Unit 


Value 




Q 


i.oo ( oOOoitU.UUUUUzu 




rljij, U i 


z4oooDy.ov5 1 oiitu.uuuuyo 


J 14 


Q 


U.U / ODDlCU.UUU^Z 


T\2 = TzA 


J 

a 


U.Ui04UltU.UUU04 


AF — /?2 /■d2 
u 






n ni Q1 fi+n nnn9n 




u.DDomu.uio 


% 





QQ fll _i_n oi 




km s ^ 


u.oy iomu-uuoo 


7 


km s ^ 


— io.iioy it u.uUi4 


c 




n 






u.uuuo 






i .U4UlllU.Uol 


-ft* 




U.yoiitU.Uif5 


log 3* 


(cgs) 


A A QO-Un Al O 


P* 


PQ 


1.211±0.U5U 




K 


5959±134 


[Fe/H] 


dex 


-0.26 ±0.10 


Mp 


Mjup 


2.303±0.068 


Rp 




1.281±0.029 


PP 


PJ 


1.096±0.067 


log gp 


(cgs) 


3.507±0.018 


a/ R, 




5.977±0.082 


a 


AU 


0.02643±0.00026 


Tp,A=o 


K 


1724±43 


Gy 


2.5tl,l 



= 6.9599 X 10** m; Mq = 1.9892 X 10^° kg; Rj^p = 7.1492 X lO'^ m; 

examine in detail the potential effects on the system pa- 
rameters of using only a single light curve. 

For survey photometry with low SNR, the durations of 
ingress and egress are ill-defined, leading to considerable 
uncertainty in the transit impact parameter and hence 
to large uncertainties in the stellar density and planetary 
radius. So-called 'follow-up' transit light curves are gen- 
erally included in the analysis of new ground-based tran- 
siting planet discoveries, and are of significantly higher 
photometric precision than the light curves produced by 
survey instruments such as WASP. Such follow-up light 
curves are typically the result of observations with a 1-2- 
m class telescope, and are of critical importance to mea- 
suring precisely basic system parameters. 

Any light curve may suffer from correlated noise, such 
as from observational systematics or from astrophysical 
sources such as stellar activity. To assess the levels of 
correlated noise in our follow-up light curves, we plot 
(Figure S]) the RMS of the binned residuals to the fit of 
each light curve as a function of bin width, along with the 
white-noise expectation. For six of our light curves, the 
RMS of the binned residuals follows closely the white- 
noise expectation, indicating that little or no correlated 
noise is present in the data. Light curves ii and vii show 
deviation from the white-noise model, however, suggest- 
ing the presence of noise correlated on timescales of ^ 1 
and ~ 10 minutes, respectively. We suggest that this red 
noise may be due to the high airmass of the target at the 
start of these observations. 



4.1. Method 
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Fig. 1. — Photometry. Upper panel: WASP-36 b discovery light 
curve folded on the orbital period of P = 1.5373653 d. For display 
purposes, points with an error greater than three times the median 
uncertainty are not shown. Lower panel: High-precision transit 
photometry, over-plotted with our best-fitting model (solid lines). 
Each individual dataset is offset in flux for clarity, and is labelled 
with a numeral corresponding to that in the first column of Table 
(2] The residuals to the best-fitting models are labelled and offset 
in the same way, and are shown in the lower part of the panel. 
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Fig. 2. — Radial velocities. Upper panel: Phase-folded radial- 
velocity measurements (Table [T] The centre-of-mass velocity, 7 = 
-13.2169 km s~l, has been subtracted. The best-fitting MCMC 
solution is over-plotted as a solid line. Middle panel: Residuals 
from the radial-velocity fit as a function of time. Lower panel: 
Bisector span measurements as a function of radial velocity. The 
uncertainties in the bisectors are taken to be twice the uncertainty 
in the radial velocities. 
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Fig. 3. — Modified Hertzsprung-Russell diagram. WASP-36 
is plotted alo ngside isochrones from the evolutionary models of 
IMarigo et aLI (2008). The isochrones are for Z = 0.01044 ([Fe/H] 
= —0.26). The leftmost curve corresponds to an age of 0.3 Gy, and 
the adjacent curve to 1.0 Gy. Thereafter the curves are regularly 
spaced at intervals of 1 Gy. 
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Fig. 4. — Correlated noise in follow-up light curves - RMS of 
binned residuals versus bin width, for light curves i - viii (solid 
lines). The white- noise expectation, where the RMS decreases in 
proportion to the square root of the bin size, is indicated by the 
straight, dotted lines. 

After modelling all available data in a combined 
MCMC analysis (see Section [331) . our 'global solution', 
we also ran several MCMCs each with just a single follow- 
up light curve in addition to the radial velocities and 
WASP photometry. Additionally, we re-ran each of these 
MCMCs applying a Gaussian prior to the stellar radius 
to impose a density typical of a main-sequence star (the 
'main-sequence constraint'). Such a constraint is usu- 
ally applied when analysing a new planet which has poor 



quality follow-up photometry (such as a single light curve 
which covers only part of the transit), and there is no 
evidence that the star is evolved or otherwise non-main- 
sequence in nature. We also performed analyses where 
the only photometry included was a single follow-up light 
curve, i.e. the WASP photometry was excluded from the 
analysis. The purpose of this is to determine whether 
the measured depth of transit is biased by inclusion of 
the WASP photometry. For these runs only, the orbital 
period was fixed to the value determined as part of our 
global solution, since this parameter is very poorly con- 
strained by a single transit light curve and a few RVs. 
The epoch of mid-transit was treated as normal, and al- 
lowed to float freely. Finally, we performed an analysis 
exlcuding all follow-up photometry; the only photometry 
analysed was the WASP data. 

4.2. Results 

We produced correlation plots between several param- 
eters, but choose to present here only plots showing im- 
pact parameter against planet radius and stellar radius 
versus stellar mass (Figures [5] and [51 respectively) . Such 
plots, whilst representative of the ensemble correlation 
plots, are particularly instructive since b and Rp are two 
of the major quantities we wish to measure, are largely 
constrained by follow-up light curves rather than by sur- 
vey photometry or by radial velocities, and can be signif- 
icantly correlated with each other, indicating a strongly 
degenerate solution. The stellar density is measured di- 
rectly from the transit light curve, and the stellar mass 
and radius, whilst interesting in themselves, are key in 
determining the values of several other system parame- 
ters of interest. 

Several conclusions can be drawn from study of Figures 
[5]and[6l and similar plots, namely: 

(1) each analysis including only a single follow-up light 
curve gives results that are consistent with our global so- 
lution, albeit with larger uncertainties. To measure the 
dispersion in the best-fitting parameter values obtained 
from each single follow-up light curve analysis, we cal- 
culated the weighted standard deviation. The standard 
deviations of b, Rp, i?* and are 0.05, 0.08 Rjup, 0.05 
Rq and 0.006 Mq, respectively. 

(2) The largest uncertainties are obtained for follow-up 
light curves that cover the smallest fraction of the transit 
(light curves i and v), as expected. 

(3) The analyses which exclude the WASP photometry 
give larger uncertainties, but these are only significantly 
so when the follow-up photometry is poor. This indi- 
cates that the WASP photometry only makes a signifi- 
cant contribution to constraining the shape and depth of 
the transit when the follow-up light curve is incomplete. 

(4) Even a partial transit light curve improves the pre- 
cision of the measured system parameters enormously 
compared to those derived solely from the WASP pho- 
tometry and the RVs. 

(5) The imposition of a main-sequence constraint does 
not significantly alter the parameters or uncertainties for 
high-precision light curves that are complete, thus indi- 
cating that WASP-36 is a main-sequence star. When the 
follow-up light curve does not well constrain the range of 
possible models, however, limiting the star to the main- 
sequence can significantly reduce the large degeneracy in 
the possible solutions. This is best illustrated by light 
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Fig. 5. — Analysis of follow-up light curves I. The MCMC posterior probability distributions for Rp and b for each of the follow-up light 
curves. The numbering of each panel corresponds to the light curve numbering in Table [2] and the l-cr and 2-cr contours are shown. In 
each case red corresponds to analysis of a single follow-up light curve plus the WASP photometry, black to the single light curve plus the 
WASP photometry with the main-sequence constraint imposed, and blue to that of a single light curve with no WASP photometry. The 
green contours indicate our global solution, and the grey contours the solution without follow-up photometry 

, and are therefore identical in each panel. 
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Fig. 6. — Analysis of follow-up light curves II. The MCMC posterior probability distributions for M* and Rt for each of the follow-up 
light curves. The numbering of each panel corresponds to the light curve numbering in Table[2]and the 1-a and 2-cr contours are shown. In 
each case red corresponds to analysis of a single follow-up light curve plus the WASP photometry, black to the single light curve plus the 
WASP photometry with the main-sequence constraint imposed, and blue to that of a single light curve with no WASP photometry. The 
green contours indicate our global solution, and the grey contours the solution without follow-up photometry, and are therefore identical 
in each panel. Also in each panel are dashed lines which are contours of constant stellar density, corresponding, from top to bottom, to 
0.7, 1.0, 1.5, and 3.0 times solar density. 
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curve (iv) , where the effects of the constraint are to de- 
crease the stellar density we find and confine the solution 
to a smaller area of parameter-space, close to the global 
solution, while largely resolving the degeneracy between 
b and Rp. 

In summary, if only one of the follow-up light curves 
had been available, we would have reached a solution 
compatible with the current best-fitting model, although 
the uncertainties on the model parameters may have 
been much greater, if the light curve was not of the 
highest precision. Obtaining additional light curves is 
clearly of benefit if one only has a light curve that par- 
tially covers transit. It is also useful to have multiple 
high-precision light curves for systems where stellar ac- 
tivity may bias the observed transit depth by varying 
amounts at difl[erent epochs, as may be the case for 
WASP-lOb ([C hristian ct al. 200^: 'Johnson ct al.' Hool; 
IDittmann et al.. 2010: Maciejcwski ct al. 2011b„ a). 

5. DISCUSSION AND CONCLUSION 

WASP-36 is a metal-poor, Solar-mass star which is 
host to a transiting planet in a 1.54 d orbit. We find 
the planet to have a mass of 2.30 Mjup, and a radius 
1.28 Rjup, meaning it is slightly denser than Jupiter. 
There is an observed corr elation between pla netary ra- 
dius and insolation (e.g. lEnoch et al.ll20lll ). with the 
more bloated planets generally receiving a greater flux 
from their star. WASP -36b is somewhat la r ger th an pre- 
dicted by the models of lBodenheimer et al.l ( 20031 ) , which 
predict radii between 1.08 Rjup(for a planet with a core 
at 1500 K) and 1.20 Rjup(for a core- less planet at 2000 

The close orbit and large radius of the planet make 
it a good target for measuring the planetary thermal 
emission, via infra-red secondary eclipse (occultation) 
measurements with, for example, Spitzer. The expected 
signal-to-noise ratios of the occultations in Spitzer chan- 
nels 1 (3.6 ^m) and 2 (4.5 ^m) are around 10 and 9 
respectively, assuming a planet with zero albedo and uni- 
form heat redistribution. 

One of the striking properties of the WASP-36 is the 
low stellar metalhcity ([Fe/H] = -0.26 ± 0.10). Giant 
planets are kno wn to be rare around s uch low-metallicity 
stars (e.g. iSantos et all l2004bl: iFischer fc Valentil 
l2005f ). although several other low-metallicity systems 
are known, including the transitin g systern s WASP-21 
([Fe/H] = -0.46 ± 0.11. .Bouchv etIdMfnoh . WASP-37 



([Fe/H] = -0.40 ± 0.12. ISimDson et aII[20Tll) and HAT- 
P-12 ([Fe/H] = -0.29 ± 0.05. iHartman erai1[200l . 

Such systems will be critical in probing our under- 
standing of the planet-metallicity correlation; proposed 
explanations for the correlation include insufficient mate- 
rial for proto-planetary cores to attain the critical mass 
needed for runaway accretion, and the suggestion that 
the high density of molecular hydrogen i n the inner galac- 
tic disk is responsible for the effect (jHavwoodl 1^009) . 
WASP-36b may also play a key role in determining 
whether stellar metalhcity is the key parameter influenc- 
ing whether or not a hot Jupiter's atmosphere exhibits 
a thermal inversio n. Insolation was in itially propounded 
as this parameter ([Fortnev et al.|[2008l) , but several plan- 
ets now appear to contradict this theory. XO-lb, for 
instance, has a relatively low insolation, and was there- 
fore predicted to lack an inversion, but Mac halek et al.l 
Hooa) report the presence of an inversion : TrE S-3b 
does not exhibit an inversion ()Fressin et ahl [2010l ) de- 
spite a pred iction to the contrar y. More recently stel- 
lar activity (iKnutson et al.ll2010[ ) and stellar metalhcity 
(jWheatlevet al.ll20lH) have been advanced as alterna- 
tives to insolation; work aiming to resolve this issue is 
ongoing. 
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